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In general correlated models, in addition to the usual adiabatic component with a spectral index 
riadi there is another adiabatic component with a spectral index nad2 generated by entropy pertur- 
bation during inflation. We extend the analysis of a correlated mixture of adiabatic and isocurvature 
CMB fluctuations of the wmap group, who set the two adiabatic spectral indices equal. Allowing 
fiadi and nad2 to vary independently we find that the wmap data favor models where the two adia- 
batic components have opposite spectral tilts. Using the wmap data only, the 2cr upper bound for 
the isocurvature fraction fiso of the initial power spectrum at kg = 0.05 Mpc~^ increases somewhat, 
e.g., from 0.76 of nad2 ~ Tiadi models to 0.84 with a prior niso < 1.84 for the isocurvature spectral 
index. We also comment on a possible degeneration between the correlation component and the 
optical depth r. Moreover, the measured low quadrupole in the TT angular power could be achieved 
by a strong negative correlation, but then one needs a large r to fit the TE spectrum. 

PACS numbers: 98.70.Vc, 98.80.Cq 



Introduction. The first studies of mixed initial con- 
ditions for density perturbations in the Ught of measured 
cosmic microwave background (CMB) angular power as- 
sumed the adiabatic and isocurvature components to be 
uncorrelated About the same time it was pointed 

out that inflation with more than one scalar field may 
lead to a correlation between the adiabatic and isocurva- 
ture perturbations 4]. If the trajectory in the field space 
is curved during inflation, the entropy perturbation gen- 
erates an adiabatic perturbation that is fully correlated 
with the entropy perturbation 0, IE S 13 • In addition, 
there is also the usual adiabatic perturbation created, 
e.g., by inflaton fluctuations. Thus, in the final angular 
power spectrum, one could have four different compo- 
nents: (1) the usual independent adiabatic component, 
(2) a second adiabatic component generated by the en- 
tropy perturbation during inflation, (3) an isocurvature 
component, and (4) a correlation between the second adi- 
abatic and the isocurvature component. In this Letter we 
assume power laws for the initial power spectra of these 
components and denote their spectral indices by rtadi, 
'^ad2j '^iso) and rtcor, respectively. Only three of these are 
free parameters, since, e.g., rtcor — (?^ad2 + ?^iso)/2. 

Although pure isocurvature models have been ruled 
out M after the clear detection of the second acoustic 
peak [lO| , the correlated mixture of adiabatic and isocur- 
vature fluctuations still remains as an interesting possi- 
bility. In 0, 0| angular power spectra have been calcu- 
lated for correlated models and compared to the CMB 
data, but the spectral indices have either been fixed or 
set equal, riadi = nad2 = ^T-iso = JT-cor- This is not neces- 
sarily well motivated theoretically. For example, if the 
entropy field is slightly massive during inflation, then 
nadi ^ 1-0 < riiso in most models. 

Recently, the Wilkinson Microwave Anisotropy Probe 
(wmap) accurately measured the temperature anisotropy 
spectrum up to the second acoustic peak [l^ and also the 



TE cross-correlation ,13], which plays an important role 
in constraining cosmological models. The wmap group 
considered the possibility of mixed models in [l3 | where, 
in order to simplify the analysis, they set the two adi- 
abatic spectral indices equal, nad2 = ^^adi- They found 
that a correlated mixture of adiabatic and isocurvature 
fluctuations does not improve the fit to the data. 

However, we would rather expect the second adia- 
batic spectral index to be close to the isocurvature one, 
'T-ad2 ~ n-iso, since both of these fluctuation components 
have been generated by the entropy perturbation dur- 
ing inflation. In this Letter we study a correlated mix- 
ture of the adiabatic and cold dark matter isocurvature 
fluctuations relaxing the "WMAP condition" by letting 
^ad2 ^adi- We show that the data clearly allows this 
and, e.g., the upper bound for the isocurvature fraction, 
/iso, slightly weakens. In this preliminary analysis we use 
the WMAP data set only, but allow riadi, ''T-ad2, and niso 
and the amplitudes of different components to vary in- 
dependently to see if there are any interesting effects. A 
more thorough analysis including other CMB and large 
scale structure data will be presented in |0| . 

Dealing with correlation. In the following we 
consider linear perturbation theory in the spatially flat 
(f2 = 1) universe. For simplicity, we assume inflation 
with two scalar fields only. Following we denote the 
inflaton field by cj) and the other scalar field by x- 

The transformation of the comoving curvature pertur- 
bation TZ and the entropy perturbation S from the time 
of the Hubble length exit during inflation to the begin- 
ning of radiation dominated era is of the form j2| 

Urad(fc)y ^0 Tssik)) \S4k)J ' ^ ^ 

where the transfer functions Tns(k) and Tss{k) carry 
all the information about the evolution of the pertur- 
bations. They are obtained by solving numerically the 
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equations of motion for the adiabatic and entropy per- 
turbations during inflation and relieating. Almost all the 
way from the generation of classical perturbations during 
inflation to the beginning of radiation dominated era the 
cosmologically interesting perturbation modes are super- 
Hubble, k <C aH. Then the evolution of perturbations is 
practically k independent . However, reheating process 
may change this conclusion. Here we assume that T^s 
is only weakly k dependent. The easiest way to deter- 
mine T-jis is to use the equation TZ = 2H9Ss/&, where a 
is the background homogeneous adiabatic field, 6s is the 
entropy field perturbation [Ss — —{2Vs/3&'^)S], and 9 is 
the angle between the inflaton field and the trajectory 
{(j){t) , xit)) in the field space as defined in 0. Writing 
'7^2 rad for the entropy generated curvature perturbation 
TTzs{k)S^{k) we get 



7^2rad(fc) = / dS^e{t)6s{k,t) 

Jt.{k) 



(2) 



where is the time when mode k becomes super- 

Hubble. If the fields and the angle 6 are evolving slowly 
during the generation of perturbations, then the main k 
dependence of TZ2 rad comes from the initial k dependence 
of (5s(fc,t*), which is very close to that of iSrad, assuming 
Tss to be also weakly k dependent. In this situation 
the k dependence of the integration lower bound 
translates into a small correction to the k dependence of 
^2 rad whcn Compared to S. This suggests that 7?.2rad 
and S have nearly identical k dependence. On the other 
hand, if the entropy field 6s is massive, it leads to iSrad 
that increases as a function of k while, at the same time 
due to slow roll, the usual adiabatic component, TZi lad = 
7?.* generated by infiaton fiuctuations, can be nearly scale 
independent as in the simplest inflationary scenarios. 

We define the correlation Cxy{k) between two pertur- 
bation quantities x and y, which in our case are TZ for 
the adiabatic and S for the isocurvature fluctuation, at 
the beginning of radiation dominated era by the formula 

{x{k)y*{k'))l^^ = ^ C,(fc) (fc - k') . (3) 
The angular power spectrum induced by the C^y will be 

C.y^^j'^C,^y[k)gZ'''''\k)9fr''\k)^ (4) 

where gi is the transfer function that describes how an 
initial perturbation evolves to a presently observable tem- 
perature (T) or polarization {E- or i?-mode) signal at the 
multipole /. 

Assuming that everything changes slowly in time and 
is weakly k dependent in a sense described after the equa- 
tion , the end result of is well approximated by the 
power laws 



(|^)'7^rad = A(|;)'"a.(fc) -f A,(A)-a.(^), 



(5) 



where and hs are Gaussian random variables obeying 

(aO=0, (a,)=0, {dr{k)al{k'))=6rs5'^''\k-k'). 

Ar, As, and B are the amplitudes of the usual adiabatic, 
the entropy generated second adiabatic, and the isocur- 
vature component, respectively. We define k — k/k^, 
where fcg = 0-05 Mpc~^ is the wave number of a ref- 
erence scale. As explained, we expect n2 — in the 
case where HO /a is nearly constant and much less than 
one during the generation of perturbations. If this condi- 
tion is not satisfied, then n2 and can take completely 
different values. General expressions for the spectral in- 
dices in terms of the slow roll parameters are derived in 
(lit. Actually, even the power law spectra (jSJ may be 
bad approximations. E.g., in double inflation numerical 
studies show that the perturbations can be strongly scale 
dependent 

Inserting (jSJ into the autocorrelations become 



Cnn = AlP^' 



and Css = B^P"^ , (6) 



while the cross-correlation between the adiabatic and 
isocurvature fluctuations is 



Cns{k)^Csn{k)^AsB~k 



"3 +"2 



(7) 



Substituting ^ and Q into and noting that the 
present total angular power is 

Ci = Cnn I + Css i + C-jis i + Csn i , 
we get for the temperature angular power spectrum 



c 



dk\ 



[Al^gl.fk^-^ + Al{gl,fk^-^ 
+ B\gl,rk'-^+2AsBgl^gl,k-^+-^ 
and for the TE cross-correlation spectrum 



(8) 



dk\ 

2„T E I2n: 



^l9Tii9n.ik^"^ + ^IdTiidnik'^"^ 
B'g'si9^ik''''+AsB {glig^^+glig'ii) fc"^+- 



(9) 



Above we defined ni, 712, and 713 so that for the scale free 
case they are zeros. To match the historical convention, 
we define new spectral indices as follows: riadi 1 = 2ni, 
riiso - 1 = 2n2, and nad2 - 1 = 2713. 

The amplitudes are not yet in a convenient form in © 
and O . The overall adiabatic amplitude at the reference 
scale fco is A^ = A'^ + A'^. Using this, the adiabatic initial 
power spectrum can be written as 



C 



A' 



(1 - y^)fc""^i-Vr''fc' 



(10) 



where = Aj/A"^ so that < < 1, Following Q 
we define the isocurvature fraction by /j^^ = (B/A)'^ and 
obtain 



C. 



ss 



A2f2 1,1 
^ Jiso'^ 
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The correlation amplitude is As B = A'^(B/A){As/A) = 
^^/isoSign(_B)\/F^ from Q. Without loss of generality, 
the total angular power spectrum can now be written as 

Ci^A\ sin^ACf'^i+cos^ACr^^^/^^Cr +/isoCosACr') , 

where > is the overall amplitude, < A < tt, and 
/iso > 0. Cf^i, Cf^2^ Cf°, and Cf"" are calculated by 
our modified version of CAMB for each cosmological 
model from ||SJ) or @ keeping A^ = Ag = B = 1. For 
example, Cf^"""""" is given by the first term in the integral 
© and the last term of integral Q gives Cf™'^^ . 

If n^di and nad2 are nearly equal or amplitude is 
close to zero or one, the adiabatic power spectrum in 
can well be approximated by a single power law 
Crni = -DA:""''"^, where D is the amplitude. How- 
ever, in the general case, an attempt to write the term 
in square brackets in (|10|) in terms of a single power 
law leads to a strongly scale dependent spectral index 
nad{k) — 1 = dlnCTin{k)/d\nk. The first derivative of 
this is always non-negative: 

rfnad(fc) ^ (1 - y^)y^(nadi - nad2)'fc"-^^+"-^^ 

The WMAP group observed that the combined CMB and 
other cosmological data favor a running spectral index 
with a negative first derivative. Thus one would expect 
that the data disfavor models where rtadi 7^ 'T'ad2, since 
this evidently leads to a positive first derivative of riad. 
However, the correlation power spectra C-ns and Csn 
may well balance the situation so that a more compre- 
hensive analysis is needed. 

Technical details of analysis. In this preliminary 
analysis we consider spatially flat (fi = 1) universe and 
use a coarse grid method leaving a sophisticated Monte 
Carlo analysis |3 for future work We concentrate 
on the neighborhood of the best-fit adiabatic model found 
in jl^. Naturally, this favors pure adiabatic models, but 
our primary interest is not to do full confidence level 
cartography here. Instead we study whether relaxing 
the WMAP constraint nad2 — J^adi has any interesting 
effects. Hence, we scan the following region of the pa- 
rameter space: reionization optical depth r = 0.11 - 0.19 
(step 0.02), vacuum energy density parameter = 0.69 

- 0.77 (0.02), baryon density ujb = 0.021 - 0.025 (0.001), 
cold dark matter density ujc = 0.10 - 0.18 (0.02), riadi = 
0.73 ~ 1.27 (0.03), nad2 = 0.55 - 1.84 (0.03), ruso = 0.55 

- 1.84 (0.03), /iso = 0.0 - 1.2 (0.04), cos A = -1.0 - 1.0 
(0.04). The best overall amplitude is found by maxi- 
mizing the likelihood for each model. We stress that, as 
in any similar analysis, the choice of the grid is a top-hat 
prior in itself. 

Since the likelihood code offered by the wmap group 
0,1^123 is far too slow for a grid method, we are able to 
use only the diagonal elements of the Fisher matrix when 
calculating the likelihoods C. Ignoring the off-diagonal 



terms increases the effective = — 21n£ by about 4 
from 1428 for well-fitted models, but since this effect is 
common to all models, it has only a small effect on the 
confidence level plots. However, we point out that the 
results presented here are mostly qualitative in nature. 

Results. From the likelihoods on the (nad2, ^^-adi) 
plane in Fig. 1(a), marginalized by integrating over all 
the other parameters, we see that the data do not espe- 
cially favor nad2 — f^adi- Clearly most of the 2a allowed 
models are in the regions where one of the adiabatic spec- 
tral indices is larger than 1.0, the other being less than 
1.0. Hence the wmap data favor models where the adia- 
batic components have the opposite spectral tilts. Using 
the full Fisher matrix of wmap our best-fit model gives 
= 1427.8 while the best-fit nad2 = '^adi model has 
= 1428.0. For comparison, our best-fit pure adia- 
batic model has — 1429.0. So, allowing for a cor- 
related mixture improves the fit slightly. However, in 
pure adiabatic models the number of degrees of freedom 
is I' = 1342 while in correlated mixed models we have 
four additional parameters leading to = 1338. Thus 
the goodness-of-fit of pure adiabatic models is about the 
same as that of mixed models: jv = 1.065 for adiabatic 
and x^ jv — 1.067 for mixed models. 

Figure 1 (b) shows that the isocurvature spectral index 
TT-iso is not limited from above. To get a constraint one 
would need to include some large scale structure data, 
which we expect to give about niso ^1-8 [l^ motivating 
our prior niso < 1.84. When the isocurvature fraction 
/iso (at fco = 0.05 Mpc^^, corresponding the multipole 
Zoff ~ 700) is large, the data favor large riiso, i.e., pos- 
itively tilted isocurvature spectrum in order to get less 
power at the smallest multipoles I. We show also by 
dashed lines how the wmap restriction nad2 = ^T-adi modi- 
fies the contours. The difference is clear in Icr region but 
2(7 regions are nearly identical. From one-dimensional, 
slightly non-Gaussian, marginalized likelihood function 
of /iso we find a 2cr upper bound for the isocurvature 
fraction, /iso < 0.84. With the restriction nad2 = "-adi, 
the bound would be about /iso < 0.76. 

Figure 1(c) should be compared with the results ob- 
tained for an uncorrelated mixture of adiabatic and 
isocurvature fluctuations in 0. Although qualitatively 
similar, the la and I0 regions of uncorrelated models 
are much smaller than in [2| due to improved accuracy 
of the data. We also notice that allowing for correlated 
models significantly enlarges the 2a allowed region in the 
parameter space. 

Since the correlation amplitude is /iso cos A, it is natu- 
ral that for a high isocurvature fraction /iso the data pre- 
fer smaller |cosA|, which is evident in Fig. 1(d). Com- 
paring one-dimensional marginalized likelihoods for riadi 
in the pure adiabatic case and in the correlated models we 
find that allowing for a correlation does not affect much 
the usual adiabatic spectrum, which is nearly scale free: 
nadi = 0.97 ± 0.06 (pure adiabatic), Uadi = 0.98 ± 0.07 
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FIG. 1: The 68.3%/lcr (white), 95.4%/2(t (hght gray), 99.7%/3<7 (medium gray), and more than 3(7 (dark gray) confidence 
levels for our general models. The best-fit model (r, Qa, (^t, ^c, J^iadi, ^^ad2, fiiso, /iso, cosA) = (0.13, 0.73, 0.025, 0.12, 1.03, 
0.64, 1.12, 0.52, -0.08) is marked by an asterisk (*) and the best-fit 'Rad2 = f^iadi model by a circle (o). The dashed lines in (b) 
are confidence levels for nad2 = riadi models and in (c) they indicate la and 2a regions for uncorrelated models, i.e., cos A = 0. 



(correlated models). However, Fig. 1(e) shows that in- 
creasing the isocurvature fraction leads to slightly stricter 
bounds for riadi. 

To demonstrate the role of the different components of 
the spectrum, we plot an angular power spectrum of a 
2(T allowed model in Fig. 2. In this particular model, a 
high 64% contribution of isocurvature to the total Cj'^ 



at the quadrupole (1=2) is allowed since the negative cor- 
relation mostly cancels the excess power. In the TE 
power spectrum there is 103% of isocurvature at the 
quadrupole. In Cj^ the cancellation between the cor- 
relation and isocurvature is not exact at the quadrupole, 
so that the isociirvature adds some power there compared 
to pure adiabatic models. The wmap measured a very 
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high quadrupole in TE spectrum, which was regarded 
as a hint for an unexpectedly large optical depth due to 
reionization t = 0.17± 0.04 [l^. In positively correlated 
models one would easily get a high TE quadrupole even 
with much smaller r. However, positive correlation is not 
favored by TT spectrum where the measured quadrupole 
is very low. In any case, in most models allowing for 
correlation adds some power at TE quadrupole, which 
explains our observation that correlated models seem to 
favor slightly smaller r than pure adiabatic models. On 
the other hand, strong negative correlation along with 
quite a large r might give a reasonable fit to both the 
low TT and the high TE quadrupole. 

The isocurvature modes introduce also another degen- 
eracy for main cosmological parameters. Namely, allow- 
ing for general initial conditions prevents one from deter- 
mining LUb from the CMB Nearly any value for ujf, 
is allowed, since it is determined by the relative heights 
of the acoustic peaks, which are also affected by even 
a small isocurvature or correlation contribution. In our 
case the degeneracy is even more severe than in , since 
we allow for independently varying spectral indices. The 
big bang nucleosynthesis calculations are valuable to de- 
termine LOb- 
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